Astronomy & Astrophysics manuscript no. cassataastroph © ESO 2008 

April 7, 2008 



5h . 
Ph. 



< 



Letter to the Editor 

GMASS Ultradeep Spectroscopy of Galaxies at z ~ 2. Ill: The 
emergence of the color bimodality at z ~ 2 

p. Cassata\ A. Cimatti^, J. Kurk^, G. Rodighiero"*, L. Pozzetti^, M. Bolzonella^, E. Daddi^, M. Mignoli^, S. Berta^, M. 
Dickinson^, A. Franceschini^, C. Halliday^, A. Renzini''', P. Rosati^' and G. Zamorani^ 

bo" 

Q ' ' Laboratoire d'Astrophysique de Marseille, OAMP, UMR6110, CNRS-Universite de Provence Aix-Marseille I, BP8, F-13376 

f~^ I Marseille Cedex 12, France e-mail: paolo . cassataOoainp . £r 

r\l , ^ Dipartimento di Astronomia, Universita di Bologna, via Ranzani 1, 1-40127, Bologna, Italy 

Max-Planck-Institut fiir Astronomie, Konigstul 17, D-691 17, Heidelberg 

Dipartimento di Astronomia, Universita di Padova, Vicolo dell'Osservatorio 2, 1-35122, Padova 
INAF-Osservatorio Astronomico di Bologna, Via Ranzani 1, 1-40127, Bologna 
CEA-Saclay, DSM/DAPNIA/Service d'Astrophysique, F-91191 Gif-sur-Yvette Cedex, France 
I — ^ ■ ^ Max Planck Institut fiir Extraterrestrische Physik, Postfach 1312, 85741 Garching bei Miinchen, Germany 

' NOAO-Tucson, 950 North Cherry Avenue, Tucson, AZ 85719, USA 
INAF-Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, 1-50125, Firenze, Italy 
f^ , '" INAF-Osservatorio Astronomico di Padova, Vicolo deirOsservatorio 5, 1-35122, Padova, Italy 

f~^ I' " European Southern Observatory, Karl-Schwarzschild-Strasse 2, D-85748, Garching bei Miinchen, Germany 

O ' Received ; accepted 

C/2 , ABSTRACT 

,^, 

Aims. The aim of this work is to study the evolution of the rest frame color distribution of galaxies with the redshift, in particular in 

the critical interval 1.4 < z < 3. 

Methods. We combine ultradeep spectroscopy from the GMASS project (Galaxy Mass Assembly ultradeep Spectroscopic Surx'ey) 

with GOODS multi-band photometry (from optical to mid-infrared) to study a sample of 1021 galaxies up to m(4.5/jm)=23. 

Results. We find that the distribution of galaxies in the (U - B) color vs stellar mass plane is bimodal up to at least redshift z = 2. 

We define a mass complete sample of galaxies residing on the red-sequence, selecting objects with log(M/Mo) > 10.1, and we study 

their morphological and spectro-photometric properties. We show that the contribution to this sample of early-type galaxies, defined 

as galaxies with a spheroidal morphology and no star formation, decreases from 60-70% at z <0.5 down to ~50% at redshift z = 2. 

^^ ' At z > 2 we still find red galaxies in the mass complete sample, even if the bimodality is not seen any more. About 25% of these red 

("^ ' galaxies at z > 2 are passively evolving, with the bulk of their stars formed at redshift z > 3. 
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1. Introduction Franzetti et al. 2007; Scarlata et al. 2007). Here we use data 

%J{ ' about spectral properties, spectral energy distributions (SED), 

d • Galaxies in the local universe show a bimodal color distnbu- morphologies and mid-IR emission that we collected for the 

■ ■ tion (Strateva et al. 2001, Hogg et al. 2002, Blanton et al. 2003). q^ASS sample to explore the content of the red sequence as 

Even If color bimodality has been observed and studied also ^ ^^^^^-^^ ^^ ^^^ ^^^^^^^ Throughout the paper, magnitudes are 

at higher redshift (Bell et al. 2004, Weiner et al. 2005, up to -^ ^j^^ ^g ^^^^ ^^^ ^^ ^^ ^^^^70 1^^^ ^-i ^ ,^ ^^_^q 3 

z ~ 1; Franzetti et al. 2006, Cirasuolo et al. 2007, up to z ~ 1.5; j q -n 7 
Giallongo et al. 2005, up to z ~ 2), no study until now has com- 
bined spectroscopic coverage and morphological analysis at z > 

1. The color bimodality suggests a different mechanism of evolu- 2. The Sample 
tion for galaxies lying on the two sequences (Menci et al. 2005; 

Scarlata et al. 2007; De Lucia et al. 2007): hence, it is exti-emely GMASS ("Galaxy Mass Assembly ultra-deep Spectroscopic 

interesting to evaluate the epoch when the color bimodality was Survey") is a project based on an ESO VLT Large Program. For 

built up. ^n exhaustive description of the survey see Kurk et al. (2008a, in 

Many authors in the literature use red sequence galax- preparation), Cimatti et al. (2008) and Halliday et al. (2008). 
ies to constrain the evolution of the early-type population The sample is extracted in the 4.5jum public image obtained 
(Bell et al. 2004a; Faber et al. 2007). Although the red peak with Spitzer Space Telescope + IRAC down to a limiting mag- 
is known to be dominated by early-type galaxies with old "itude of OT4.5 =23.0, and contains 1021 galaxies in a field of 
passive stellar populations, a contamination of star forming 6.8'x6.8'. The selection at A.5^im is more sensitive to stellar 
galaxies with colors reddened by dust, or by eai'ly-type spi- mass and less affected by dust extinction than optical bands, 
rals is also present (Bell et al. 2004b; Cassata et al. 2007; Moreover this selection produces a "negative" k-correction at 
z > 1 .4, as the peak of the stellar spectral energy distribution 
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Fig. 1. Rest-frame (L'^ - B) color versus stellar mass in six red- 
shift intervals. A separate panel is dedicated to the forming clus- 
ter at (z) - 1.61. The diagonal continuous line indicates fit to 
the red sequence, as the dotted ones show the scatter around the 
best fit, for the different redshift bins. The dashed line indicates 
10.1 M/Mq, the mass limit at which the red sequence is com- 
plete at any redshift in our sample. Color-coded symbols indi- 
cate galaxies in different morphological classes: red, blue, green 
and cyan symbols represent respectively early-types, spirals, ir- 
regulars and undetected objects. In each panel the histogram of 
colors is also reported. The bottom-right panel reports the color- 
mass relation for all the galaxies, regardless of the redshift. 



With the aim of studying the redshift range 1.5 < z < 3, the 
epoch when the crucial processes of massive galaxy formation 
took place, ultra deep spectroscopy was carried out using FORS2 
on galaxies pre-selected with a cut in photometric redshift of 
Zphot >l-4. The spectroscopic coverage was complemented with 
available literature redshifts, reaching a completeness of about 
50%. For the rest of the population, accurate photometric red- 
shifts are available (Kurk et al. 2008a). The median of the galaxy 
redshift distribution is (z) = 1.2, with 190 objects with spec- 
troscopic redshift z > 1.4, mostly coming from GMASS spec- 
troscopy. There are several spikes in the redshift distribution, the 



most dominant being at (z) = 1 .61, discussed in a parallel paper 
(Kurk et al. 2008b). 

The photometric SEDs were derived for all the galax- 
ies in the sample using the public images available in the 
GOODS-South ai-ea in 11 bands: optical (ACS-HST, BVIz, 
Giavalisco et al. 2004), near-infrared (ESQ VLT-hISAAC, 
JHKs) and mid-infrared (Spitzer+IRAC, 3.6/zm, 4.5/zm, 5.6/zm 
and 8//m). We used the synthetic spectra of Maraston (2005; 
M05) to fit these SEDs, adopting a Kroupa IMF, limiting the 
fit to Aresi < 2.5/im. We use exponentially declining star for- 
mation histories, with timescales spanning from 0. 1 to 30 Gyr, 
plus the case of constant star formation rate. Models with ages 
between 0. 1 Gyr and the age of the Universe at the redshift of 
each galaxy are retained in the best-fit procedure. Extinction is 
a free parameter in the optimization, using the extinction curve 
of Calzetti et al. (2000). The fitting procedure minimizes i\\t x^^ 
and the best fit model gives an estimate for the age, the e-folding 
time of the SFR r, the extinction Ay and the stellar mass. At 
the same time, absolute magnitudes in Johnson UBV bands are 
derived. The dataset is complemented with the Spitzer-MlPS 
data pubUcy available for the GOODS-South/GMASS region, 
to check for possible activity signs in the 24yL/m data. 

The high resolution imaging provided by ACS-HST allowed 
an accurate visual classification of all the galaxies in the sample. 
The analysis has been performed independently by two of us (PC 
and GR) on the ACS band closest to the rest-frame B-band. The 
classification scheme is based on 4 classes: 1. spheroidal galax- 
ies (ellipticals, SO and compact objects); 2. spirals; 3. irregular 
galaxies; 4. undetected in the optical bands (and thus not classifi- 
able). On the basis of this analysis, we classified 198 spheroidal 
galaxies, 496 spirals and 269 irregulars, while 58 objects are un- 
detected in the optical bands. 

At redshift z > 1 .2, where the ACS z-band maps the blue UV 
light in the rest-frame, surface brightness dimming and morpho- 
logical k-correction effects may be important. However, the sep- 
aration between early- and late-types should remain robust, as - 
though fainter in the UV- ellipticals remain symmetrical. Spirals 
instead can appear more late type, as red bulges get fainter, the 
surface brightness of disks dimms, but knots of star formation 
brighten in the UV. 

3. Color bimodality 

Precise redshift measurements are extremely important to reduce 
uncertanties in the SED fitting procedure, to derive robust abso- 
lute magnitudes, stellar masses, and other SED parameters. The 
GMASS spectroscopic coverage in the range 1.4 < z < 2.5 al- 
lows on one hand a proper calibration of the photometric red- 
shifts, and, on the other hand, provides us with high-quality fits 
of the spectral energy distribution for each galaxy. The opti- 
mal measure of the rest-frame magnitudes allows us to check 
for the bimodality in colors in a redshift interval that has not 
been much explored until now. In Figure [1] we report the rest- 
frame iU-B) color versus stellar mass for different redshift bins, 
as well as for the global population. Each redshift bin contains 
roughly the same number of objects. In each panel a small inset 
shows the color distribution at that redshift. The forming cluster 
at (z) = 1.61 is excluded from the bins at z > 1.2, but we ded- 
icated a separate panel to it. The quality of our sample allows 
us to claim that the color distribution, looking at both the U - B 
histograms and the (U - B)vs M/Mq plane, is bimodal at least 
up to z = 2, as well as for the global population. We note that 
galaxies show also a sort of bimodality in mass, with red galax- 
ies being more massive of the blue ones. This is partly due to 
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a selection effect: at a given redshift, the completeness limit is 
higher for red than blue galaxies. Conversely, on the other side 
of the diagram, the lack of massive blue galaxies is a real effect. 

To characterize and define the red sequence, first we se- 
lect objects redder than the valley in the global color dis- 
tribution, namely (U-B)=l.l. Then, we fit the global color- 
mass distribution of these galaxies with a straight line (U- 
B)=fl-h/7(log(M/M0), obtaining for the slope a value /7=0.0943. 
Finally, we use this slope to fit the red sequence in each redshift 
bin, leaving only the intercept a as a free parameter The result of 
this procedure is overplotted on the data in Fig.[T] Interestingly, 
the intercept a does not evolve significantly with redshift: we 
measure a A(f/ - B)--Q.2 between z ~0.5 and z ~2.5. This is 
consistent with the evolution expected for a galaxy formed in a 
instantaneous burst at redshift z- 5. 

In Figure [T] we color coded galaxies of different morpholo- 
gies. We find that 70% of the spheroidal galaxies reside on the 
red sequence. The remaining 30% show blue colors. On aver- 
age, red spheroidal galaxies are more massive than blue ones 
(the median of the mass distribution is 2 x 10'" and 5 x 10^ M© 
respectively for red and blue). On the other hand, 77% of late 
type galaxies have blue colors, with the remaining 23% lying on 
the red sequence. At redshifts z > 2 many objects have colors 
z-K> 2, being very faint or undetected up to the ACS z-band. 
For these no morphological analysis could be made. 
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Fig. 2. Comparison between Specific Star Formation Rates 
(SSFR) coming from SED fitting and from 24jum luminosity, 
for galaxies at z < 2. Red points represent galaxies lying on the 
red sequence, having log(M/M0) > 10.1 (the mass completeness 
limit at any redshift), while blue empty circles show galaxies in 
the blue cloud. Empty triangles are galaxies with passive spec- 
tra, empty diamonds are galaxies with passive spectra but some 
signs of emission lines. Downward arrows identify galaxies for 
which the 24/im luminosity is just an upper limit, as the they 
have S 24^m < 25yuJy. 



4. Red sequence composition 

In general, galaxies can have a red color because they have old 
stellar populations, or because their star-formation activity is ob- 
scured by dust that reddens their colors. Hence, the red sequence 
is well known to be a mix of galaxies with different proper- 
ties (Bell et al. 2004b; Cassata et al. 2007; Scarlata et al. 2007; 
Franzetti et al. 2007). Here we have the opportunity to explore 
the morphological, photometric and spectroscopic properties of 
red galaxies up to redshift z ~ 3. 

To separate the blue and red population in each redshift bin, 
we use the best fit to the red sequence and add an offset of 
Am = -0.15 magnitudes to end up on the valley between red 
and blue galaxies. This offset is comparable to the scatter of the 
objects around the fit to the red sequence. To allow a compar- 
ison between galaxies at different redshifts, we define a mass 
complete sample, selecting only galaxies on the red-sequence 
having logM/M© >10.1. This is roughly the smallest mass for 
red galaxies at redshift z ~ 3 in our sample. With this criterion, 
we select a sample of 197 galaxies. The aim of this section is to 
study the properties of these galaxies, by combining morpholog- 
ical information and SED analysis, including also 24yum data. 

In Figure|2]we study the star formation properties of galaxies 
in this sample, comparing the Specific Star Formation Rates (the 
star formation rate for unit mass, SSFR) estimates coming from 
SED fitting and from 24/vm photometry. We plot only galaxies 
up to z - 2, because beyond this limit the data show to much 
scatter. We report also blue galaxies for comparison. The total 
infrared luminosity (8-lOOOyum) was computed by rescaling var- 
ious star forming templates to the 24/im flux of each source (see 
Daddi et al. 2007 for details). Then, this total luminosity was 
converted into a Star Formation Rate using the relation given 
by Kennicutt (1998). We use 25fj.iy as a conservative detection 
limit for 24jum sources in GOODS (corresponding roughly to 
S /N < 2), reporting as upper limits the values of the SSFR be- 
low this limit. 

The agreement between the two SSFR estimates up to z = 2 
is quite good, in particular for blue galaxies, even if a large scat- 



ter is present. The two measurements still correlate for red mas- 
sive galaxies, considering also that the 24jum SSFR provides an 
upper limit for the bulk of them. The large majority of the mass 
complete red sample has log(SSFRs££)) < -1.2 and no reliable 
detection at 24jum. The quite good correlation between the two 
SSFR and the agreement with spectroscopic properties ensure 
that this diagram can be used to identify galaxies which have 
very low levels of SSFR and are likely to be the ones with the 
lowest or no star formation. In particular, we define passive those 
galaxies having log(SSFRs£fl) < -1.2: at this rate, a galaxy 
would take about 15 Gyr to double its mass. To be more con- 
servative, we excluded galaxies having a reddening Ay >1 and 
those with log(SSFR24;„„) > -1 and S24fjm > 25/iJy (that are re- 
liable 24/im sources). We checked the position on the diagram 
of the 16 galaxies in the mass complete sample having passive 
GMASS spectra (see Cimatti et al. 2008 for the 13 at z > 1.4): 
14 of them occupy the "passive" region of the diagram, while 
2, that also show some emission lines in their spectra, are in the 
"active" region. 

At z > 2 the scatter about the relation becomes significant 
and we therefore decided not to plot galaxies with z > 2 on 
Fig. |2j this scatter is probably due to increased uncertanties on 
the 24//m luminosity and/or to an increasing number of dusty 
star forming objects. In this range we decided to use the same 
limits to define passive galaxies, which exclude objects with 
large 24//m luminosities. Finally, among the 197 galaxies be- 
longing to the mass complete sample, 103 (~ 50%) could be 
considered passive according to these criteria. 

It is interesting to combine the results of the morphologi- 
cal classification with those of the analysis of the star formation 
activity, for our mass complete sample of galaxies on the red se- 
quence. We find that all the objects classified as spheroidals are 
actually passive galaxies. This means that at we do not find, at 
any redshift, red sequence massive ellipticals with a secondary 
episode of star formation. On the other hand, not all the galax- 
ies that are passive according to the criteria mentioned above 
are morphologically spheroidals. In fact, about 20% of the 103 
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passive galaxies on the red sequence are morphologically late- 
types. The large majority of them are spirals dominated by a 
bright bulge. 

We summarize these results in Fig. [3] where we show 
the contribution to the massive part of the red sequence 
(log(M/M0 > 10.1)of galaxies that show little or no signs of star 
formation (log(SSFR5£o) < -1.2&Av < 1 & log(SSFR24;im) < 
- 1) and have an early-type morphology. Since at z > 2 we do not 
have a morphological classification for all the galaxies, at that 
redshift we combine SED and lAfiva information to define galax- 
ies with low SSFR. In the literature, authors either include or ex- 
clude bulge dominated spirals as a part of the early-type family, 
we therefore show here both cases. It can be noted that the trend 
of the two cases is similar, with the contribution of passive early- 
type galaxies mildly (or not at all) decreasing with the redshift 
from z ~0 to z ~ 2. The fluctuations in the frequency could be 
due to cosmic variance: i.e, the lack of massive galaxies on the 
red sequence in the bin around z ~ 1 may imply that our survey 
samples an underdense region of the Universe at that redshift 
(see the photometric redshift histogram in Kurk et al. 2008a). 
The results at redshift z < 1 .5 are in good agreement with pre- 
vious results: Franzetti et al. (2007) find an increasing contri- 
bution of star forming galaxies in the red sequence between 
z ~0 and z ~ 1.2. Moreover, Scarlata et al. (2007) found that 
the fraction of morphologically early-type galaxies (including 
bulge dominated spirals) in a sample of photometrically selected 
ETGs decreases from -60% at z = 0.3 to ~ 45% at z = 0.9 . 
Renzini (2006) reported that at z~ 58% of the red galaxies are 
morphologically early-type, similar to our findings. At z~ 0.7 
Bell et al. (2004b) and Cassata et al. (2007) find respectively 
75% and 66% of early-types contributing to the red sequence: 
these results can be easily conciliated with ours considering that 
the former include Sa galaxies, while the latter consider pure el- 
lipticals/SO. 

Here we reach higher redshift. We show that even at 1 .5 < 
Z < 2 the color distribution is bimodal, with the massive part 
of the red sequence dominated by passive spheroidal galaxies 
(Fig.lT^O. Even though at z > 2 bimodality in the color distri- 
bution is no longer evident, many galaxies do show red colors 
and masses above log(M/M0) - 10. 1. We find that ~ 25% of 
them have low or no signs of star formation (compared to 60% at 
z < 2), having log(SSFR5£o) < -1.2 & log(SSFR2v„) < -1 & 
Ay < 1. Half of these red and massive galaxies have a spheroidal 
morphology, even if the remaining half can not be morpholog- 
ically classified, because they are barely detected or undetected 
in the ACS z-band. Nevertheless, all of them have BzK colors 
compatible with being passive (Daddi et al. 2004). Combining 
their ages (greater than 0.5 Gyr) to their redshift, we can con- 
clude that this population of massive galaxies must have formed 
the bulk of their stars during a burst at redshift z > 3. It is plau- 
sible to argue that the bimodality at z > 2 vanishes as a result of 
the decreasing number of passive red galaxies and of the conse- 
quent increase of dusty star forming objects. 

Finally, we emphasize that, even at z < 1 and more signifi- 
cantly at z > 1, the fraction of massive late-type galaxies lying 
on the sequence is significant, and evolves with redshift. Thus, 
it could be misleading to identify red-sequence with passive old 
galaxies. This must be taken into account when red sequence 
galaxies are used to constrain the evolution of the luminosity 
and mass function for the passive early-type population, as in 
Befl et al. (2004a), Scarlata et al. (2007) and Faber et al. (2007). 




Fig. 3. Contribution to the massive part of the red sequence 
(log(M/Mo) > 10.1) as a function of the redshift, for early-type 
galaxies with little or no signs of star formation. Red strip in- 
cludes bulge dominated spirals, while blue dashed one includes 
pure spheroidal galaxies. 



5. Conclusions 

We study the rest-frame (U - B) color distribution for a sample 
of 1021 galaxies selected at 4.5/im. We took advantage of the 
ACS/HST high resolution imaging and multiwavelength cover- 
age spanning from ACS bands to Spitzer IRAC mid-infrared and 
lAjim to characterize morphologies, colors and SED properties 
of galaxies in the sample. 

1 . We find that the classical color bimodality is preserved up to 
at least z = 2. This can be seen both looking at the galaxies 
in the (U-B) vs Mass plane, and at the {U - B) distribution. 
The presence in the sample of a forming cluster at (z) - 1.61 
does not enhance the bimodality at 1.5 < z < 2. 

2. The massive part of red sequence (log(M/M0) > 10.1) is a 
mix of early- and late-type galaxies at all redshifts. The con- 
tribution of early-type galaxies to the red sequence is about 
60-70% at z ~ 0.5 and about 50% at z = 2, depending on the 
inclusion or exclusion of bulge dominated passive spirals in 
the early-type class. 

3. Even if the color bimodality vanishes at redshift z > 2, still 
there are red galaxies with log(M/Mo) > 10.1. About 25% 
of these are passively evolving, with ages >0.5 Gyr. This im- 
plies that we find a rich population of massive galaxies with 
redshift of formation z > 3. 
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